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Abstract —In this work, we propose an abductive framework 
for biosignal interpretation, based on the concept of Temporal 
Abstraction Patterns. A temporal abstraction pattern defines 
an abstraction relation between an observation hypothesis and 
a set of observations constltntlng its evidence support. New 
observations are generated abductlvely from any subset of the 
evidence of a pattern, building an abstraction hierarchy of 
observations in which higher levels contain those observations 
with greater interpretative value of the physiological processes 
underlying a given signal. Non-monotonic reasoning techniques 
have been applied to this model in order to find the best 
interpretation of a set of initial observations, permitting even 
to correct these observations by removing, adding or modifying 
them in order to make them consistent with the available domain 
knowledge. Some preliminary experiments have been conducted 
to apply this framework to a well known and bounded problem: 
the QRS detection on ECG signals. The objective is not to 
provide a new better QRS detector, but to test the validity of an 
abductive paradigm. These experiments show that a knowledge 
base comprising just a few very simple rhythm abstraction 
patterns can enhance the results of a state of the art algorithm 
by significantly improving its detection Fl-score, besides proving 
the ability of the abductive framework to correct both sensitivity 
and specificity failures. 

I. Introduction 

A prominent objective of biosignal processing research is 
to provide classification algorithms for identifying the inter¬ 
esting phenomena from signal samples. Such classifiers can 
be grouped into two main approaches: (1) a knowledge-based 
approach, which aims to model the domain or, alternatively, 
to model an expert reasoning process m; and a (2) learning- 
based approach, which builds a model by estimating the 
underlying mechanisms that produce the data of a training 
set m. Once the classifier is obtained and validated, it behaves 
at a logical level as a deductive system from a data vector 
representing the signal. The present work departs from the in¬ 
trinsic limitations of the deductive framework for coping with 
the interpretation of biosignals. Indeed, deductive approaches 
apply a monotonic consequence relation, so that any conclu¬ 
sion cannot be retracted as new evidence becomes available. 
This entails a propagation of errors from the first processing 
stages onwards, narrowing the capability of making a proper 
identification as new processing stages are successively added. 
Usually, both the above-mentioned approaches overcome this 
weakness through an artificial adoption of a casuistry-based 
strategy, yielding to unsatisfactory results so far. 


The present work starts from the hypothesis that temporal 
abductive reasoning El provides a more appropriate frame¬ 
work for the computational interpretation of biosignals. Its 
greatest strength lies in its non-monotonic nature, so that 
the conclusions are guessed as conjectures inferred from the 
available evidence at each time, and further information, at 
different levels of abstraction, can modify these conclusions. 
For example, a deductive arrhythmia classifier cannot correctly 
identify an ECG fragment in which the beats were not properly 
detected. In contrast, an abductive interpreter could conjecture 
the presence and morphology of a beat from its context, like 
a human can reconstruct a speech despite failing to identify 
all its constituent sounds. 

This paper provides an abductive framework for biosignal 
interpretation on the basis of the notion of Temporal Abstrac¬ 
tion Pattern, which is a knowledge representation formalism 
that defines an abstraction relation among a set of domain 
observables, structuring the problem domain in an abstraction 
hierarchy, where higher abstraction levels convey more seman¬ 
tic content and more interpretative results. Patterns are gen¬ 
erated by means of attributed regular grammars El, extended 
in order to address the representation of temporal knowledge 
through the Simple Temporal Problem formalism Q, previ¬ 
ously adopted for temporal abduction in diagnosis B- 

The rest of this paper is outlined as follows; Section [I^ 
defines the main notions of the proposal. Section III describes 
an interpretation algorithm, based on heuristic searching. Sec¬ 
tion IV shows the representation of basic ECG knowledge, and 


how it is used for QRS detection. Finally, section [V] discusses 
validation results with respect to a state of the art algorithm, 
and section [VT] provides some conclusions. 


II. Definitions 

A Simple Temporal Problem (STP) El defines a set of 
temporal constraints between pairs of temporal variables, Ti 
and Tj. A temporal constraint is a closed interval L{Ti, Tj) = 
[dijjbij] restricting the admissible values for the difference 
Tj — Ti, as Uij < Tj — Ti < bij. Formally, a STP can be rep¬ 
resented as a tuple N = (T, C), where T = {Ti, ..., T„} is a 
set of temporal variables, and C = {L{Ti,Tj);l < i,j < n} is 
a set of temporal constraints between them. A tuple (ti,..., f„) 
is called a solution of N if the assignment Ti = fi,...,T„ = 
satisfies all the constraints. 




The basic representation entity of our abductive framework 
is the observable, which is formally defined as a tuple 
q = {rj, A, Tb, Te), where 77 is the name of the observable, A = 
{Ai,An ) is a set of attributes to be valued, and Tf, and Te 
are two temporal variables representing the beginning and the 
end of the observable. We denote by Q = {qq, qi,qn} the 
set of observables of a particular domain, comprising the vo¬ 
cabulary to describe the phenomena of interest. An observable 
can be observed in multiple instances called observations. An 
observation is defined as a tuple o = ( 77 , u, 4 , fg), where 77 is 
the name of the observable being instantiated, 17 is a value 
assignment for each attribute of the observable, and tb and te 
are the specific values for the temporal variables Tb and Te- 
We denote by 0{q) = {of,..., o®,...} the set of observations 
of the observable q. An example of an observable is g = 
(QRS, shape, Tft, Te), which represents a QRS complex with 
a single attribute describing its morphology. This observable 
may be instanced in different observations, like for example 
of = (QRS, shape = QS, 00:32.123,00:32.201). 

Observables are related through abstraction relations defined 
by temporal abstraction patterns. A temporal abstraction pat¬ 
tern models the necessary knowledge to allow the conjecture 
of an observation of a high-level observable from observations 
of lower level observables. We provide a procedure for dy¬ 
namically generating abstraction patterns, based on the formal 
language theory. The set Q of observables can be considered as 
an alphabet represented by their corresponding names. Given 
an alphabet Q, a formal grammar G denotes a pattern of 
symbols of the alphabet, describing a language L{G) C Q*, 
as a subset of the set of possible strings of symbols of the 
alphabet. Let G°'P be the class of formal attributed grammars 
of abstraction patterns. A grammar G G G“p is syntactically 
defined as a tuple (Vat, Vt, H, R). The production rules in R 
are of one of the following forms: 


H — qu —> g[^]G 

G ^ q[l]D I q[l] I A 

H = qn is, the initial symbol of the grammar, and it plays 
the role of the hypothesis guessed by the pattern. Vn is the 
set of non-terminal symbols of the grammar. Vr is the set of 
terminal symbols of the grammar, gathering together: a set of 
observables Qg ^ Q. being q ^ qn for all q G Qg, that can 
be abstracted by the hypothesis; a set of temporal descriptions 
[/] in the form of conjunctions of constraints between the 
temporal variables of the observable produced by the rule and 
all the observables previously generated; and the empty string, 
that is represented by A. Some simple examples of grammars 
to describe different basic cardiac rhythms in terms of QRS 


complexes are detailed in section IV 

Given a grammar G G G“p, we provide a constructive 
method for representing a set of abstraction patterns Pq = 
{Pi,..., Pi,...}. An abstraction pattern shows a temporal ar¬ 
rangement between a set of observables, possibly appearing 
repeatedly, to be abstracted by a new observable qn- We call 
findings to these occurrences of observables that lead to qn. 


which informally behave as observations that have not yet been 
observed, that is, predictions generated by the grammar defin¬ 
ing an abstraction pattern. Thus, Mp = {m{, m^,..., = 

{rj, A,T^,T ^),...} is the set of findings of the observable q in 
P. Pg gathers together the set of abstraction patterns that share 
the same observable qn to be abstracted, so they represent the 
different ways of hypothesizing qn- 

Definition 1: Given G G G“p, a temporal abstraction 
pattern P = {qnp, Mp, Np,Ilp) consists of a hypothesis 
qpp', a set of findings Mp = U^gCg 
evidence supporting qp ^; a temporal network Np between the 
temporal variables involved in qpp and Mp, defined by the [1] 
attributes of G; and an observation procedure lip to compute 
the attribute values of qpp from the observed evidence. 

Below we show how an abstraction pattern is built by 
following the productions of a grammar. In every step, a new 
observable is added, and a set of temporal constraints among 
this finding and those generated before is introduced in Np. 

1) Symbol H entails initializing an abstraction pattern: 

P ^ {qpp ,Mp = 0, {Tp = {Tb",Tf}, Lp = {L(Tf, rf)})) 

2) All those productions H — qpp q[l]G entail: 
P^{qHp,Mp = {ml},{TpU{T^,T^},CpUJZ{P,ml))), 

C{P,ml) = L{Tl,Tl) U {L{Ti,Tj)\Ti G Tp A Pj G {tI,tI}} 

3) All those productions G — q[l]D \ q[l] entail: 

P ^ {qHp,Mp U {ml}, {Tp U {T^, re"}, Cp U C{P, ml))), 

C{P,ml) = L{T^,T^) U {L{T,T^)\T gTp A G [T^,T^}} 

In case a temporal descriptor is omitted in a production of the 
G°-P grammar, it is assumed the ’after’ relationship between 
the new manifestation and the set of previous manifestations. 

For instance, all those productions G —> qD \ q entail: 

P^{qH,MpVJ {ml}, {Tp U {T^, Te"}, £p U C{P, ml))), 

C{P,ml) = L{TtT^)yj{L{T,Tt) C Z+\T G {Tl,Ti}Amj G Mp} 

Definition 2: Let Q be a set of observables and V a set 
of abstraction patterns. We say V induces an abstraction 
relation in Q x Q, denoted by qi\< qj if and only if there 
exists a pattern P G V such that: 

1) qi G Mp 

2 ) qj = qHp 

3) qi^^qi, where |<+ is the transitive closure of |< 

The abstraction relation allows us to say that “qj abstracts 
qi”. This gives a hierarchy structure to the observables of a 
domain, allowing us to define an abstraction model from a set 
of abstraction patterns. Fig. illustrates this hierarchy with 
the observables of the knowledge base defined in section 
Definition 3: We define an abstraction model as a tuple 
Ml = {V, Q, |<) where V induces an abstraction relation |< 
over a set of domain observables Q. 

Definition 4: We define an interpretation problem as a 
tuple IP = {0,M), where O = {oi,..., o^,...} is a set of 



observations requiring interpretation, and A4 is an abstraction 
model of the domain. 

An interpretation problem so defined gives a different status 
to the evidence with respect to abductive diagnosis, where 
the explicit difference between normal and faulty behaviors 
leads to the definition of faulty findings Q. Only when a 
faulty finding is provided, the diagnostic process is triggered. 
In contrast, an interpretation problem based on an abstraction 
model gives all the findings the same status. The objective is to 
provide a new description in terms of those observables with 
the highest possible abstraction level. To do so, and under 
certain constraints, an observation is tried to be assigned to 
a specific finding of an abstraction pattern, so that a new 
observation is obtained for the hypothesis of the pattern. 

Definition 5: Given an interpretation problem IP, a match¬ 
ing relation for a pattern P S 7^ is an injective relation in 
Mp X O, defined by ^ Oj, iff m| = (p, A, T^, TJ) € Mp 
and Oj = G O. 

From the notion of matching relation we can design a 
mechanism for abductively interpreting a subset of observa¬ 
tions in O through the use of abstraction patterns. Thus, a 
matching relation for a given pattern allows to hypothesize new 
observations from previous ones, and to iteratively incorporate 
new evidence to the interpretation by means of a hypothesis 
generation-and-test cycle. The notion of abstraction hypothesis 
defines those conditions that a subset of observations must 
satisfy in order to be abstracted by a new observation. 

Definition 6: Given an interpretation problem IP, we de¬ 
fine an abstraction hypothesis as a tuple h = {oj^, Ph,^h), 
where Ph G V, •^h'^ ^Ph ^ ^^d we denote Oh = 

{oji ,..., Ojn} = image{^h), satisfying: 

1 ) o^ G 0{qHp) 

2 ) ^IVp{Oh) 

3) satisfy the constraints of Np. 

An abstraction hypothesis assigns a set of observations to 
the findings of the pattern, giving them the role of evidence for 
the hypothesis. Even though the matching relation is a matter 
of choice, and therefore a conjecture by itself, some additional 
constraints could be assumed as default assumptions. An 
important default assumption in the abstraction of periodic 
processes states that consecutive observations are related by 
taking part of the same hypothesis, defining the basic period 
of the process. This assumption leads to consecutive findings 
of an observable within an abstraction pattern to be matched 
to consecutive observations during the abstraction task. 

As a result of an abstraction hypothesis, a new observation 
is generated, that should be included in the set of domain 
observations, so O = O L) {of^}. The abstraction process 
then iterates until no new observations are generated. The 
set of observations that may be abstracted in an interpretation 
problem IP is 0{domain{^)), that is, the set of observations 
corresponding to observables playing the role of findings in 
some abstraction pattern. Figure [T] shows an example of an 
abstraction hypothesis of the Extrasystole pattern. 

Definition 7: Given an interpretation problem IP, we de- 



Fig. 1. Example of an abstraction hypothesis of the Extrasystole pattern. 
At the bottom of the figure the five initial beat annotations are shown. Four 
of these observations are abstracted by QRS observations that are matched 
with the evidence of the Extrasystole pattern, allowing to discard the 
incorrect observation. The top of the figure displays the constraint 

network of the pattern. 

fine an interpretation as a temporally consistent set of ab¬ 
straction hypothesis I = {/ii, hm}, satisfying OOh- = 

0 , i ^ j- 

The restriction Oh^ H Oh^ — 0 establishes that in the 
same interpretation an observation cannot be abstracted by 
more than one hypothesis, formalizing the notion of alternative 
hypotheses as those abstracting the same observation. An inter¬ 
pretation can be rewritten as / = {Of, Pi,<—i, Nj), where: 
of = {of, is the set of conjectured observations; 

Pj = {Pi,Pm} is the set of abstraction patterns used in the 
interpretation; ^/=^i U...U {^Pi U U Mp^) x O 

is the global matching relation; and Nj = Ni U ... U Nm. is 
a network between all the temporal variables involved. The 
temporal assignment made by is a solution of V/. 

III. Interpretation algorithm 

The purpose of the interpretation algorithm is to obtain the 
best interpretation for a given set of initial observations. Since 
a number of different sets of abstraction hypotheses may be 
obtained with the same base evidence, we need a method to 
discriminate some interpretations against other. To evaluate 
and compare the generated interpretations, we rely on the 
following general principles: The coverage principle, which 
prefers interpretations explaining more initial observations, 
and that is calculated as C{I) = | 0;i|/|0(domain(|< 

))|; the simplicity principle, also known as Occam’s razor, 
which prefers interpretations with fewer hypotheses, and that 
is calculated as S{I) = 1/(1 -f \Of\)', the abstraction 
principle, which prefers interpretations using terms of higher 
abstraction levels; and the predictive principle, which prefers 
interpretations that properly predict future observations. 

Once we are able to evaluate and compare interpretations, 
solving an interpretation problem can be posed as a heuristic 
search on the space of consistent interpretations. The proposed 
search strategy uses coverage as the main heuristic function, 
and considers simplicity when the coverage is equal in more 
than one interpretation. The main drawback of coverage as 
heuristic is that is non-admissible, and therefore optimality can 
not be guaranteed, so we require an algorithm efficient with 













Algorithm 1 Partial Expansion K-Best First Search algorithm. 
1: function PE-KBFS(7P, A) 

2: var/° = (0,0,0,0) 

3: SET_FOCUS(/°, miTlti (Oi € O)) 

4: var open = sorted([({1.0, 1.0), 7°)]) 

5: var closed = sorted([]) 

6: while open 7^ 0 do 

7: for all I € open[0 ... 77] do 

8: r = NEXT(GET_SUCCESS0RS(7)) 

9: If not 7' then 

10: closed = closed U {({1 - C(7), 1/S(7)), 7)} 

11: else if C{I') = 1.0 then 

12: return 7' 

13: end if 

14: var val = (1 - C[I'),l/S{r)) 

15: open = open yj {{val, I'), {val, I)} 

16: end for 

17: end while 

18: return min{closed) 

19: end function 


this type of heuristics and that saves memory and computing 
time. We propose an evolution of the K-BFS method 0, 
named Partial Expansion K-BFS, whose pseudocode is shown 
in Algorithm [T] The algorithm takes as input an interpretation 
problem IP, and returns the first found interpretation with 
full coverage, or the interpretation with the highest coverage 
factor and highest simplicity if no full coverages are found. 
The algorithm also needs a K parameter, which determines 
its exploratory capabilities, and which is set to the maximum 
number of abstractions that can be made to an observable: 
K = max{{\{qj \ gi|< qj}\,qi G Q})- The intuition behind 
this decision is that at a certain point in the interpretation, and 
with equal coverage factor, we give an opportunity to every 
possible hypothesis to continue the interpretation. 

The distinctive features of PE-KBFS with respect to KBPS 
are that (1) in each node expansion, only one successor is 
obtained, greatly limiting the growth of the open list; and 
(2) the heuristic estimation of the newly generated nodes 
is propagated to the parent node, resulting in a significant 
reduction in the number of expansions provided that child 
nodes are generated ordered by their valuation. These two 
features are aimed at reducing the number of node expansions, 
which in our case results in fewer generated interpretations. 

The search process starts with the so-called trivial interpre¬ 
tation 7°, the one with no abstraction hypotheses. At each step, 
each of the K best interpretations generated so far is evolved 
in a hypothesis-and-test cycle. This cycle is driven by the 
concept of focus of attention, which is established according 
to the predictive and abstraction principles and enables to 
combine different inference modes based on the state of the 
interpretation. The focus is initially set to the earliest initial 
observation in the interpretation problem, and its evolution 
depends on the state of the hypothesis-and-test cycle, managed 
by the GET_SUCCESSORS() function, detailed in Algorithm]^ 

In the hypothesis step, the focus is on an observation not 
yet abstracted. Then, abductive reasoning is applied and a 
hypothesis of a higher abstraction level is generated. This is 


Algorithm 2 Method for obtaining the successors of an 
interpretation in the adopted reasoning scheme. 

1: function get_successors(7) 

2: var focus = get_eocus(7) 

3: var succ = 0 

4: if is_observation(/ocus) then 

5: succ = abduce(7, focus) 

6: else if is_finding(/ocus) then 

7: succ = subsume(7, focus) U deduce(7, focus) 

8 : end if 

9: return succ 

10: end function 

11: function subsume(7,/) 

12: var succ = 0 

13: for all Oi € O I / <— Oi do 

14: I' = {Of,Pi,^iU f^Oi,Ni) 

15: succ = succ U {7'} 

16: end for 

17: return succ 

18: end function 

19: function deduce(7, /) 

20: var succ = 0 

21: for all [P G V \ f e 0{qHp)} do 

22 : h = {o"^,P,<D) 

23: 7' = (Of U {of^}, P, u {P}, U / ^ of^, Af 

24: 0 = OU{of^} 

25: for all £ next_findings(P) do 

26: I" = {0’^.,Pv,^v,Nv) 

27: SET_FOCUS(7",m‘') 

28: SUCC = succ U {7"} 

29: end for 

30: end for 

31: return succ 

32: end function 

33: function abduce(7, Oi) 

34: var succ = 0 

35: var V°' = [P GV \ OiG 0{qi), qi G Mp} 

36: for all P G V°' do 

37: for all nrp' G do 

38: h = {o^’^ ,P,mS' Oi) 

39: 7' = (Of U{of^},P/U{P},^/ 

40: 0 = OU{of^} 

41: SET_FOCUS(7',of^) 

42: succ = succ U {7'} 

43: for all G next_findings(P) do 

44: I" = {Of,,Pp,^p,Np) 

45: SET_FOCUS(7",m‘^) 

46: succ = succ U {7"} 

47: end for 

48: end for 

49: end for 

50: return succ 

51: end function 


implemented in the ABDUCEO function, which explores all the 
patterns that may abstract the focused observation and, for 
each of them, it transfers the focus to the evidence findings 
predicted by the pattern, switching to the test step of the cycle. 
This reasoning mode enforces the abstraction principle, since 
it always generates a hypothesis of high abstraction level on 
the basis of an observation of a lower level. 

In the test step, the focus is on a finding predicted by 










an abstraction pattern of the current interpretation. In this 
case, two possible inference modes may be used, subsumption 
and deduction. The SUBSUME() function looks for an existing 
observation that can be consistently matched with the focus 
finding, while the DEDUCE() function creates a new observation 
hypothesis for the matching and then looks for the lower level 
evidence supporting it. Since an unmatched finding can be seen 
as a prediction of the knowledge involved in an interpretation, 
these two inference modes enforce the predictive principle, but 
inasmuch as subsumption leads to simpler interpretations it is 
preferred over deduction. 

The NEXT_EINDINGS() function, used both in abduction and 
deduction, allows to extend an existing hypothesis, returning 
the possible next findings that may generate a pattern P, 
obtained through the analysis of the Gp grammar. If this 
function returns nothing, then the focus of attention is set 
to the earliest observation in the domain of the abstraction 
relation with the highest possible abstraction level, starting a 
new hypothesis-and-test cycle. 


IV. An application to ECG QRS detection 


In order to test the practical feasibility of this framework, 
we have defined a simple knowledge base to tackle a well 
known and bounded problem: QRS detection on ECG signals. 
This knowledge base assumes as initial observations the QRS 
annotations identified by any other algorithm, and describes 
some of the most common rhythm patterns in terms of 
temporal distances between beats. The temporal consistency 
with any of these patterns allows us to discriminate those 
annotations that do not really correspond with QRS complexes, 
as well as to proactively search for those QRS that, being not 
initially annotated, are predicted by any of the patterns. 

The knowledge base defines the following set of observables 

Q = { 90 , 91 , 92 , 93 , 94 , 95 }, with: 


qo = (BeatAnn, 0, Tb = Tf, = T) qi = (QRS, 0, Tb = Te = T) 
q 2 = (NormalRhythm, 0, T{,, Te) 53 = (Bradycardia, 0, T;,, Te) 
g 4 = (Tachycardia, 0, Ti,, Te) qs = (Extrasystole, 0, T;,, Te) 


qo represents a beat annotation, the starting point of the 
interpretation process, and that is considered an instantaneous 
observable, hence = T^', qi represents the QRS complex, 
which is also an instantaneous observable; 52 , 93 54 

are different observables for regular cardiac rhythms, and the 
only difference between them is the maximum and minimum 
admitted distance between consecutive beats; 95 represents an 
extrasystole, which is a cardiac rhythm characterized by the 
presence of a premature beat and the subsequent return to the 
previous rhythm with a compensatory pause. 

Below we show the grammars and observation procedures 
relating these observables (all temporal constraints are ex¬ 


pressed in milliseconds): 


Gpg :H = BeatAnn 
Gp, :H ^ QRS 
Gp^ :H = (N|B|T)Rhythm 
A 
B 

Gp,, :H = Extrasystole 

G 

D 

E 


—5-BeatAnn[(i] 

—)>QRSo[^ 2 o] a 
^QRSi[/2i] B 
—>QRS„[/2n] B I QRS„[Z2n] 
—>QRSq[^5o] G 
^qrSi[;5i] d 

—>QRS2[^52] E 

-AQRS3[(53] 


['ll 

[« 20 ] 

[^ 21 ] 

[hn] 


[*50] 

[^ 51 ] 

[* 52 ] 

[^ 53 ] 


{^{,^QRS_,pBeatAnn) ^ [_150,150]} 

= [ 0 , 0 ]} 

{NRR = [475,1333], BRR = [1000, 2000], TRR = [200, 600], 
^ (N|B|T)RR} 

{RR = - tWS„_2^ 

^ (N|B|T)RRn [0.5 ■ RR, 1.5 ■ RR], 

= [ 0 ^ 0 ]} 

{^(jnQRSo jnExtrasystole^ _ j-Q | 

= [ 200 , 2000 ]} 

{RRO = , 

L(T®®i,T‘’R'^2) = [200,0.9 ■ RRO]} 

{RRl = , 

L(T‘^R®i,T‘’R®3) = [1.7-RRO, 2.3-RRO], 

L(T®®TT‘’R'^3) = [1.25-RRl, 4-RRl] 

jiExtrasystolej _ j-Q 


Ho = maxt{ip‘^{ECG{t))),t & L(ro,T®*****™)} 

Hi = maxt{\ECG(t) - mode{ECG)\),t £ I,{To,r®®)} 

In the procedures, Tq represents the time origin, in this case 
the start of the ECG recording, and therefore L(To,T) is the 
interval of possible assignments for the variable T. 

The first pattern Pq is a purely deductive pattern, since 
it does not abstract any observable. It is included in the 
knowledge base in order to allow the discovery of new beat 
annotations in those cases in which a rhythm pattern predicts 
the presence of an annotation, but it is not in the initial 
evidence. The observation procedure IIo uses a wavelet-based 
method detailed in ||9l, performing a wavelet transformation ip 
of the signal in the predicted interval, and emits as hypothesis 
the time point in which the energy of the transform is greater. 

The Pi pattern abstracts the initial annotations in observa¬ 
tions of QRS complexes. The pattern introduces as temporal 
constraint [Zi] a maximum distance between the beat annota¬ 
tion point and the QRS location. The observation procedure 
Hi establishes the time instant of the QRS in the point of 
maximum deviation of the signal with respect to the baseline, 
estimated as the mode of the signal in the observation interval. 
As shown in Eig. [T] this procedure allows to change the timing 
of the initial annotations. 

Regular rhythms ( 72 , 93 and (74 share the same pattern gram¬ 
mar P 2 , which requires the presence of at least 3 consecutive 



QRS complexes with a bounded distance between them. For 
NormalRhythm, this distance is established to a heart rate 
of 45-120 bpm; for Bradycardia, the admitted rate is 30- 
60 bpm; and for Tachycardia, 100-300 bpm. The overlap 
between patterns allows the correct interpretation of those 
fragments that are on the frontier of two rhythms. In addition 
to these static constraints, [l 2 n] introduces an additional con¬ 
straint limiting the instant rhythm variation. These patterns do 
not have an observation procedure If, since the temporal limits 
are set to the time point of the first and last QRS complexes. 

The extrasystole pattern P 5 is also defined only by temporal 
constraints. An extrasystole requires the presence of exactly 4 
consecutive beats. The first two define the reference rhythm, 
while the [Z 52 ] constraint sets an upper bound on the duration 
of the second RR interval, forcing the third beat to be pre¬ 
mature. [^ 53 ] introduces the compensatory pause constraints, 
requiring the difference between the second and fourth beats to 
be approximately two times the previous RR, and a signihcant 
increase in the RR interval between the third and fourth beats. 

With this simple knowledge base, we are able to apply the 
interpretation algorithm presented in section |ng in order to 
obtain the best interpretation of a set of beat annotations, 
which may be obtained by any external algorithm. To avoid an 
excessive exploration of the search space if no interpretations 
with full coverage factor were found, a restriction in the 
search procedure was included, consisting in limiting the 
size of the open list to at most K elements (in this case, 
K = 4) whenever the interpretation process exceeds the real 
duration of the interpreted signal fragment. This ensures that 
the interpretation is performed in soft real time. As a result of 
the interpretation, a new set of annotations will be generated 
from the QRS observations present in the best interpretation 
obtained by the algorithm, provided that the observation is 
included in a rhythm pattern. 

V. Validation results 

To validate the ability of the presented framework to revise 
and correct the results of classical deductive approaches, we 
have used a state of the art algorithm providing a set of 
beat annotations for an ECG recording. These annotations are 
converted to instances of the BeatAnn observable, and form 
the initial observations upon which the abductive interpreta¬ 
tion process is applied. The selected algorithm is the WQRS 
algorithm m, included in the standard distribution of the 
WFDB software package ca. The validation dataset consists 
of a selection of ECG recordings showing regular rhythms 
and extrasystoles, and comprises all the 18 recordings of 24 
hours duration of the Normal Sinus Rhythm (NSR) database 
and 20 recordings of 30 minutes duration from the MIT-BIH 
Arrhythmia database, from the Physionet initiative ca. 

Table 1 ^ shows the comparative results of the original algo¬ 
rithm and the corrected output through abduction, in terms of 
sensitivity, positive predictivity, and the combined El-score. 
As can be noted, the abstraction process can slightly decrease 
the sensitivity, but is compensated with an improvement in 
the positive predictivity to get in the majority of cases an 



WQRS 

WQRS + Abduction 

Record 

Se 

P-H 

FI 

Se 

P-H 

FI 

NSR-16265 

100.00 

99.74 

99.87 

99.97 

99.76 

99.86 

NSR-16272 

98.24 

89.79 

93.83 

97.74 

93.66 

95.66 

NSR-16273 

99.99 

99.93 

99.96 

99.95 

99.98 

99.96 

NSR-16420 

99.98 

99.79 

99.88 

99.92 

99.92 

99.92 

NSR-16483 

99.98 

99.88 

99.93 

99.98 

99.97 

99.97 

NSR-16539 

99.97 

99.79 

99.88 

99.76 

99.92 

99.84 

NSR-16773 

99.99 

99.96 

99.97 

99.96 

99.95 

99.95 

NSR-I6786 

100.00 

99.97 

99.98 

100.00 

99.98 

99.99 

NSR-16795 

99.99 

99.87 

99.93 

99.96 

99.88 

99.92 

NSR-I7052 

99.98 

99.52 

99.75 

99.95 

99.68 

99.81 

NSR-I7453 

99.98 

99.72 

99.85 

99.94 

99.89 

99.91 

NSR-I8I77 

99.98 

99.63 

99.80 

99.91 

99.74 

99.82 

NSR-I8I84 

99.99 

99.55 

99.77 

99.98 

99.79 

99.88 

NSR-I9088 

100.00 

98.29 

99.14 

99.98 

98.37 

99.17 

NSR-I9090 

99.99 

99.70 

99.84 

99.99 

99.79 

99.89 

NSR-I9093 

100.00 

99.87 

99.93 

99.99 

99.88 

99.93 

NSR-I9I40 

100.00 

99.82 

99.91 

100.00 

99.84 

99.92 

NSR-I9830 

99.93 

98.74 

99.33 

99.86 

98.99 

99.42 

NSR-Gross 

99.90 

99.08 

99.49 

99.83 

99.40 

99.61 

MIT-100 

100.00 

99.95 

99.97 

99.95 

99.95 

99.95 

MIT-101 

99.93 

99.80 

99.86 

99.80 

100.00 

99.90 

MIT-102 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 

MIT-103 

99.94 

100.00 

99.97 

100.00 

100.00 

100.00 

MIT-104 

100.00 

97.58 

98.78 

99.25 

98.93 

99.09 

MIT-105 

99.72 

91.72 

95.55 

98.00 

97.24 

97.62 

MIT-107 

99.89 

98.13 

99.00 

99.66 

99.94 

99.80 

MIT-108 

99.59 

86.20 

92.41 

95.54 

97.79 

96.65 

MIT-109 

99.86 

100.00 

99.93 

100.00 

100.00 

100.00 

MIT-IIl 

99.94 

99.94 

99.94 

99.77 

100.00 

99.88 

MIT-112 

100.00 

99.91 

99.95 

100.00 

100.00 

100.00 

MIT-113 

100.00 

100.00 

100.00 

100.00 

99.67 

99.83 

MIT-115 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 

MIT-117 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 

MIT-122 

100.00 

99.95 

99.97 

100.00 

100.00 

100.00 

MIT-123 

100.00 

100.00 

100.00 

99.84 

99.06 

99.45 

MIT-209 

100.00 

99.53 

99.76 

100.00 

99.88 

99.94 

MIT-212 

100.00 

99.91 

99.95 

100.00 

100.00 

100.00 

MIT-230 

100.00 

99.95 

99.97 

100.00 

100.00 

100.00 

MIT-234 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 

MIT-Gross 

99.95 

98.58 

99.26 

99.62 

99.64 

99.63 


TABLE I 

QRS DETECTION CORRECTION RESULTS 


advance in the combined El-score. Still, cases such as records 
MIT-103 or MIT-109 show that the abstraction process is also 
able to correct sensitivity failures. In order to prove that the 
improvement on the detection is signihcant, the Wilcoxon 
statistical test has been applied to the differences on the El- 
score, obtaining a p-value of 0.008 on the records of the 
Normal Sinus Rhythm database, a p-value of 0.033 on the 
records of the MIT-BIH Arrhythmia database, and a combined 
p-value of 0.001 on the full set of records. This demonstrates 
that the improvement is statistically signihcant. 

VI. Conclusion 

In this paper we have presented an abductive framework 
to address the problem of biosignal interpretation. We have 
conducted a validation study to prove that the non-monotonic 
behavior is able to review the results of classical approaches, 
trying to enhance a state of the art QRS detector using a 
very simple knowledge base. The objective was not to provide 
a new algorithm for QRS detection, but to prove that an 


















































abductive approach can overcome the limitations of deductive 
methods. Results demonstrate the capacity of our proposal to 
correct and improve the outcomes of the original algorithm. 
Future plans include building a more complete knowledge 
base, modeling other ECG constituents, like P and T waves, 
explicitly considering the morphology of the significant waves, 
and including more rhythm patterns. This is expected to 
provide a more interpretative information, helping to improve 
previous results on ECG processing at all abstraction levels, 
from QRS detection to arrhythmia recognition. 
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